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ABSTRACT 

In a previous paper we have proposed a non-Newtonian, phenomenological description of the effective 
gravitational force acting between exotic dark matter and baryons, which was shown to fit weU the 
kinematics in the inner regions of low surface brightness, gas-rich galaxies. The Gravitational Suppres- 
sion (GraS) scheme is, admittedly, an ad-hoc model; it does however successfully address interesting 
cosmological issues. Here, we test the potential generality of the GraS concept. GraS predictions on 
motions on scales measurable in dwarf spheroidal galaxies and in the solar neighborhood show significant 
departures from those of Newtonian mechanics. By analyzing such systems, we find that the paradigm 
of a universal GraS potential cannot be rejected. Indeed, compared to Newtonian predictions, GraS 
provides a better description of data, when realistic dark matter density profiles are considered. 

Subject headings: Cosmology: dark matter — Galaxies: dwarf — Galaxies: halos — Galaxy: solar 
neighborhood — Gravitation — Stellar dynamics 
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1. INTRODUCTION 

Dark matter (DM) has proved to be an essential compo- 
nent of the self-consistent cosmological framework, within 
which observations on a wide range of cosmic scales are 
generally interpreted. Despite its pervasive role, however, 
its nature is still elusive and its properties mostly specu- 
lative. In an attempt to gain insight of its fundamental 
characteristics, and motivated by the apparent conflicts 
faced by Cold Dark Matter (CDM) scenarios when com- 
pared with observations on small cosmological scales (e.g., 
Sellwood & Kosowsky 2001; Ostriker & Steinhardt 2003), 
we have recently hypothesized that the gravitational at- 
traction between DM and baryons may be suppressed on 
the kpc scale (Piazza & Marinoni 2003, hereafter paper 
I). The gravitational suppression (GraS) hypothesis was 
suggested by the simple observation that all the various 
inconsistencies between theory and observations appear to 
be characterized by a common physical length. GraS min- 
imally modifies the existing paradigm, i.e. it does not alter 
large-scale interactions or baryon and DM self-dynamics. 

In paper I we modeled the effective force acting between 
visible and dark particles by adding a short-range Yukawa 
contribution to the standard Newtonian potential. We 
also investigated the consistency and universality of such 
an empirical gravitational paradigm, fixing the model pa- 
rameters using observations of the rotation curves of low 
surface brightness (LSB), gas-rich galaxies, the cores of 
which are thought to be DM-dominated. The scale length 
A over which the Yukawa contribution is effective is about 
one kpc, and its strength on small distances is equal and 
opposite to the Newtonian one. In other words, the in- 
ner dynamics of LSB galaxies seems to point towards a 
"total suppression" scenario, whereby DM is completely 
oblivious of normal matter in its immediate vicinity. 

The analysis of paper I was mainly concerned with ro- 
tationally supported disk galaxies having circular velocity 



of order ~ 100 km/s. Here, we investigate consequences 
and predictions of the proposed gravity model on a class 
of systems having different matter distributions and dy- 
namical symmetries, with characteristic velocities nearly 
a factor 10 smaller than those of typical LSB galaxies. In 
particular we aim: 

a) to assess the ability of the model in describing 
the kinematics of dwarf spheroidal (dSph) galaxies, i.e. 
in a physical regime of spherically symmetric, pressure- 
supported, low-velocity, virialized systems. This study is 
motivated by recent claims (e.g., Lokas 2002) according to 
which a universal Navarro, Frenk & White (1997, NEW) 
DM profile can reproduce the shape of the velocity dis- 
persion profiles of dSph galaxies only if a substantial (and 
poorly justified) amount of tangential anisotropy in the 
velocity distribution of dSph stars is assumed. By con- 
fronting kinematical data of dSph galaxies with the GraS 
predictions, we wish to verify whether the latter achieves 
the same interpretative success as with LSB galaxies 

b) to check the ability of GraS in describing both kine- 
matics and space distribution of stars perpendicular to 
the galactic disk in the vicinity of the Sun (e.g., Oort 
1932; Siebert, Bienayme, & Soubiran 2003). Since the 
presence of a spheroidal, non baryonic, DM halo embed- 
ding the thin galactic disk affects the density distribu- 
tion of stars near the galactic plane (Bahcall 1984), it 
is imperative to quantify the dynamical effects on stel- 
lar motions which are the consequence of GraS. It has 
been noted (Kuijken & Gilmore 1989) that an intrinsi- 
cally different non-Newtonian theory (Modified Newtonian 
Dynamic: MOND, by Milgrom 1983), while capable of ex- 
plaining the flat rotation curves of galaxies without requir- 
ing the presence of a dark matter halo, leads to inconsis- 
tent predictions on the vertical component of the gravita- 
tional acceleration experienced by solar neighborhood disk 
stars. Does GraS score any better? 
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2. THE MODEL 

As in paper I, we model the proposed short-range mod- 
ification of gravity by adding to the standard Newtonian 
potential a repulsive Yukawa contribution which is active 
only in the mixed (dark-visible) sector. We can then 
write the total gravitational potential acting on baryons 
as i> = i>jv — ^Y, where satisfies the usual Poisson 
equation 



y *w = 47rG ipB+ Pd) 



(1) 



where pB and pn are respectively the baryonic and DM 
densities. The Yukawa contribution typically dies out 
exponentially with a radial scale length A, according to 

(y^ - X-^)<i>Y ^'iirGpo- (2) 

Notice that (f>y is sourced by the DM energy density only. 
The resulting gravitational potential between a DM and a 
baryonic particle at a distance r from each other goes like 
^g-r/A _ iyj._ In what follows we assume A = 1 kpc as 
suggested by the study of LSB galaxies in paper I. 

3. MASS DISTRIBUTION AND DYNAMICS IN DSPH 

GALAXIES 

Dwarf spheroidal galaxies, being small and DM dom- 
inated, provide a test for the validity of GraS in a low- 
velocity, shallow potential regime. Since the typical lu- 
minous dimension of a dSph is ^ 1 kpc, its baryonic dy- 
namics is expected to be strongly influenced by the spe- 
cific form of our gravitational potential, which is signifi- 
cantly dampened with respect to the pure Newtonian one 
on these scales. Thus: given a realistic DM density profile, 
is the baryon dynamics derived within the GraS paradigm 
in agreement with the observations of the stellar velocity 
dispersion in dSph galaxies? 

The stellar radial velocity dispersion profile ct^ induced 
by the total potential $ can be obtained by solving the 
Jeans equation for a pressure supported spherical system 



dr 



{pb<jI) 



„/3 2 

2-pB<Jr + PB-T = 

r dr 



(3) 



where (3 = 1 — (Tg{r) /a^ (r) accounts for the anisotropy in 
the stellar velocity field. 

The luminous component of a dSph galaxy is well de- 
scribed by the Sersic profile (Sersic 1968; Ciotti 1991) 

/(r) = /oexp[-(r/ro)i/™] (4) 

with exponential index ~ 1 (but see Caldwell 1999) where 
I{r) is the surface brightness. The baryonic energy density 
PB{r) can be obtained from I{r) by deprojection: 

dl{t) dt 



PB{r) 



TT Jr 



dt Vi^ - r'^ 



(5) 



(see Mazure & Capclato (2002) for an explicit anal}d;ical 
solution of this integral) where T* is the baryonic mass-to- 
light ratio. To first approximation, T* is independent on 
the radius for the objects considered. The solution (5) can 
be applied to equation (3) and in the inversion of equa- 
tion (1) to obtain ^m- The DM component po is also 
needed, as it contributes to both the Newtonian (1) and 
the Yukawa (2) potentials. We describe po via a general- 
ized NFW profile of inner slope index 7 

= (r/r.)ni + r/r.)3-7' 



where po is the critical density of the Universe, (5 is a char- 
acteristic overdensity and a radial scale parameter. Ac- 
cording to N-body simulations, 7 = 1 (Navarro, Frenk & 
White 1997) and 7 = 1.5 (Moore et al. 1998) bracket the 
plausible range of variation of this parameter (Navarro et 
al. 2003). 

The virial radius r-v of an object is that within which the 
mean density is A ~ 200/9o ^ Q.023{M /h-^MQ^^^ kpc. 
The radius marks the transition from the outer to 
the inner r^"' density behavior. A concentration param- 
eter for the DM halo c = r^/rg can be defined, which in 
N-body simulations exhibits a mild dependence on the to- 
tal mass of the halo. We have not found explicit results for 
objects of masses < lO^^Af©, so we extrapolate the results 
of Eke, Navarro & Steinmctz (2001) to small masses and 
use, for the NFW 7=1 case, the formula 



c-Y=i = 81 



M 



h-'^Mp 



-0.07 



(7) 



while, for the steeper profile 7 = 1.5 we set c^=i 5 = 
(l/2)c^=i (Jing & Suto 2000). 
Equation (2) in the case of spherical symmetry reduces 

to 

d'^^Y 2 d<^Y 



dr2 



rdr #=4-^^-W- («) 



and in paper I an integral expression for <^y in terms of 
Pd was outlined. However, for generic 7 5^ 1 only numeri- 
cal solutions of (2) can be found, while the solution of the 
Poisson equation (1) can always be given in closed ana- 
lytic form for any value of 7. Once the total potential $ is 
obtained, equation (3) can be solved for the radial veloc- 
ity dispersion. A more useful quantity is the line-of sight 
velocity dispersion ct;os, which is related to ar through 
(Binney & Mamon 1982) 



T*7(r) 



r2\ pB{t)a^{t)t 



dt (9) 



and provides a direct comparison with observations. 



4. MODELS VS. OBSERVATIONS 

Spatially resolved kinematic data became recently avail- 
able for two dSph galaxies: Fornax (Mateo 1997) and 
Draco (Kleyna 2002). Under the assumption that these 
galaxies are in dynamical equilibrium, as found by Ma- 
teo (1998); Piatek et al. (2002); Odenkirchen et al. (2001); 
Kleyna et al. (2002)), we can compare the observed veloc- 
ity dispersions with those predicted by GraS. Assuming 
that the measurement errors of the dSph velocity profiles 
are normally distributed (possible systematics such as the 
influence of binary populations are negligible as shown for 
example by De Rijcke & Dejonghe (2002)), we compare the 
observations of aios with Equation (9) and constrain the 
scale parameters of the DM density profile via least square 
minimization. We also fix the parameters characterizing 
the stellar distribution (i.e. the baryons properties) as 
given by the observations, as shown in Table 1. 
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Fig. 1. — The observed velocity dispersion profile of Fornax (Ma- 
teo 1997) is shown together with four different best fitting velocity 
dispersion profiles inferred assuming isotropic orbits (/3 is kept fixed 
to in the fit). Theoretical predictions are computed in the case 
of a DM inner density profile slope of 1 (Navarro, Frenk & White 
1997) or 3/2 (Moore et al. 1998) and according to both the pure 
Newtonian and the GraS paradigms, as labeled. Error bars indicate 
1(T uncertainties. 

Numerical minimization is performed assuming as free 
parameters of the model the virial mass M of the dark halo 
embedding the dSph and the velocity anisotropy parame- 
ter P which, for simplicity, is assumed to be independent 
of the radial scale. Thus, our results do not rely on any a- 
priori assumed anisotropic model: the velocity anisotropy 
is self-consistently evaluated by our fitting scheme. The 
resulting best fitting values are presented in Table 2. 

Our results may be described as follow: 

a) the predicted GraS gravitational acceleration gener- 
ated by a two component fluid of baryons (distributed ac- 
cording to the Sersic profile) and DM (distributed accord- 
ing to the generalized NFW profile) is in agreement with 
the stellar velocity dispersions observed in the dSph galax- 
ies Fornax and Draco, if the stellar orbits are isotropic as 
shown in Figs. 1 and 2 

b) the best fitting velocity anisotropy parameter /3 ~ 
derived within the GraS gravitational paradigm is consis- 
tent with the results of N-body simulations of DM halos 
(e.g., Thomas et al. 1998). Table 2 shows that a Newto- 
nian analysis of the kinematics within a cuspy, NFW DM 
profile agrees with the observations only if the existence 
of a strong tangential anisotropy in the dSph stellar or- 
bits is assumed (see also Lokas (2002)). The amount of 
tangential anisotropy needed to fit the profiles is, how- 
ever, unphysically large if confronted with observations 
of stars in the solar neighborhood (e.g., Chiba & Berrs 
2000) or the evidence that elliptical galaxies are mod- 
erately, radially anisotropic (e.g., Gerhard et al. 2001) 



Fig. 2. — The same as in fig. 1 but for the Draco velocity 
dispersion profile. Data are from Kleyna et al. (2002) 



c) the best fitting mass within the virial radius depends 
on the assumed concentration parameter (eq. 7). For re- 
alistic values of c, and for all gravity models, the mass of 
Fornax and Draco is larger than commonly assumed (see 
Stoehr et al. (2002) for a similar conclusion). Since in 
pressure supported systems there is a well known degener- 
acy of density profiles versus velocity anisotropy (Binney 
& Tremaine 1987), the best fitting masses present a non 
negligible degree of covariance with respect to the best 
P estimate. The more negative the /3 parameter (as in 
the Newtonian case), the bigger is the inferred mass (see 
Table 2). Moreover, the steeper the cusp of the profile 
(higher 7), the more tangential is the velocity distribution 
required to fit the data, in agreement with the finding of 
Lokas (2002). We also note that our results for Mt/Ly 
(the mass-to-light ratio inside the dSph tidal radius rt) 
agree in order of magnitude with other independent es- 
timates based on Newtonian analysis (e.g., Kleyna et al. 
(2001,2002); Odenkirchen et al. (2001)) 

d) dSph galaxies have low densities and therefore they 
lie in the regime of small accelerations for which MOND 
provides an alternative interpretation of the observed kine- 
matics, without reliance on the DM hypothesis. However, 
Lokas (2001) and Kleyna et al. (2001) found that the best 
fitting values of Oo, the characteristic acceleration scale 
(ao ~ lO^'* cm s^^) below which the laws of Newtonian 
dynamics change, are different for different dSph (low for 
Fornax, much higher than expected in the case of Draco). 
They concluded that the Draco poses a serious problem for 
MOND. At variance, the GraS parameters are stable and 
robust even in the case of these two systems which have 
similar velocity dispersion but very different luminosities. 
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5. THE SOLAR NEIGHBORHOOD DYNAMICS 

We next test GraS predictions in the sub-kiloparsec, 
small velocity regime of stc^Uar kinematics in the direction 
z perpendicular to the galactic disk. In the vicinity of the 
disk the gravitational potential can be locally described by 
(Bahcall 1984) 



= 47rG(pg^^ + pr°)- 



(10) 



In what follows, derivations with respect to z and r are 
meant while keeping the other coordinates - cylindrical 
and spherical respectively - fixed. Moreover, equation (10) 
is derived under the assumption that the rotation curve is 
flat (/9^^'° oc r~^) which is, to a good approximation, a 
valid assumption at our position in the Galactic disk. We 
also note that in the solar neighborhood the halo density 
is about one tenth of the disk baryonic density (Gates, 
Gyuk & Turner 1995). As a consequence, the halo effects 
on the nearby stellar dynamics can be treated as a per- 
turbation, approximated by expanding in e = p^^° / p'b^^ 
(Bahcall 1984). The dynamics of a single isothermal pop- 
ulation of density p with velocity dispersion cr^ is then 
described by the first moment of the Boltzmann equation 
a1dp/dz = —pd^/dz. 

We next estimate the correction terms introduced by 
GraS into eq. (10). Note that, since is spherically 
symmetric. 



2 (f<^Y 

dz"^ r dr \r/ dr^ 



d^<^Y 1 







(11) 



where the height from the disk z is small with respect to 
our distance r from the center of the Galaxy. In order 
to estimate the first term on the right hand side of (11) 



we use equation (8). Since the distance of the Sun from 
the center of the Galaxy is about 8.5 kpc, we expand the 
solution in the r ^ X limit to obtain 



$^(r) = -47rGAV5?'°(r-) 



1 + 



(12) 



Deviations from the Newtonian equation (10) are therefore 
expected at order (A/r)^ in the halo density contribution. 
For instance, in the isothermal approximation cx r~^, 
we obtain 



92$ 



1 + e 1 



2X' 



(13) 



At the Sun's location, the correction to the Bahcall param- 
eter e is thus negligible and of order ~ 3%. We conclude 
that GraS is fully compatible and does not compromise 
our current understanding of the kinematics in the solar 
neighborhood. 

GraS offers a phenomenologically viable paradigm with 
which to reconcile mismatches between observations and 
CDM predictions on small cosmological scales. 
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Table 1 
Dwarf parameters 



Parameters 


Fornax 


Draco 


ro(kpc) 


0.35 


0.18 


n(kpc) 


2.85 


0.95 


m 


1 


1.2 


distance (kpc) 


138 


82 


Lv{Lq) 


1.55-10^ 


2.6- 10^ 


T^Mq/Lq) 


2.5 


2.5 



Note. — All data are taken from 
Mateo (1998) except the Sersic scale ra- 
dius, the exponential index m of the lu- 
minosity density profile and the tidal 
radius of Draco which have been re- 
cently rc-cstimatcd by Odcnkirchcn ct 
al. (2001) using SDSS data. The bary- 
onic mass-to-light radius represents the 
average value for a typical globular- 
cluster star population (e.g., Mateo et 
al. 1991; Pryor 1996) 



Table 2 
Best Fitting parameters 





DM inner 
slope 


Gravity Model 


P 


M 
(lO^Me) 


Mt/Lv 




Fornax 


7=1 


GraS 


„+0.2 
U-0.4 




0.65 


0.96 




7 = 3/2 


GraS 


-n 

"•J--0.2 


9.51^ 


0.65 


0.95 




7=1 


Newton 






0.23 


0.91 




7 = 3/2 


Newton 


q-l-0.7 


13.6t| 


0.76 


0.92 


Draco 


7 = 1 


GraS 


n 9-1-0.18 
"•^-0.23 




16 


0.78 




7 = 3/2 


GraS 






15 


0.80 




7 = 1 


Newton 




■■■•'-0.7 


5 


0.78 




7 = 3/2 


Newton 




00-^80 
^'^-50 


33 


1.07 



Note. — All the uncertainties represent la errors 



